Collective quantum jumps of Rydberg atoms 
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We study an open quantum system of atoms with long-range Rydberg interaction, laser driving, 
and spontaneous emission. Over time, the system occasionally jumps between a state of low Rydberg 
population and a state of high Rydberg population. The jumps are inherently collective and in fact 
exist only for a large number of atoms. We explain how entanglement and quantum measurement 
enable the jumps, which are otherwise classically forbidden. 
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Perhaps the strangest aspect of quantum mechanics is 
the notion that merely observing a system changes it. 
This concept is taken to the extreme in the quantum 
Zeno effect, where the constant observation of a system 
inhibits a transition that would otherwise take place 
Another equally striking phenomenon is quantum jumps, 
where a system under continuous monitoring occasion- 
ally switches between two distinct states [1, [!]■ Quan- 
tum jumps have been observed in many settings, such as 
trapped ions 0-0 , photons 0] , electrons [1| , and super- 
conducting qubits [9|. In these experiments, the object 
being observed is a single particle or can be described 
by a single degree of freedom. But the recent interest in 
generating multi-particle entanglement IB" 121 raises the 
question of how large systems of entangled particles be- 
have under constant observation. For example, do they 
undergo collective quantum jumps? 

In this paper, we show how entanglement and quantum 
measurement lead to collective quantum jumps of Ryd- 
berg atoms. A Rydberg atom is an atom excited to a high 
energy level n. The dipole-dipole interaction between two 
Rydberg atoms is strong and allows one to entangle many 
atoms over long distances This interaction has at- 
tracted recent interest for quantum-information process- 
ing 13- 18 1 and many-body physics 19l-24|. 

We consider a group of atoms laser-driven to the Ry- 
dberg state and spontaneously decaying back to the 
ground state. Classical mean-field theory predicts two 
stable collective states, one with low Rydberg popula- 
tion and one with high Rydberg population. Classically, 
the system should remain in one of the stable states. 
However, we find that quantum fluctations drive tran- 
sitions between the states, resulting in quantum jumps. 
The jumps are inherently collective and exist only for a 
large number of atoms. Our results may be extended 
to other settings, such as coupled optical cavities 25-28j| 
and quantum-reservoir engineering i29--32|. 

Two atoms in the same Rydberg level experience an 
energy shift V due to their dipole-dipole interaction [III ■ 
The dependence of V on inter-particle distance R can 
take several forms. In the presence of a static electric 
field, V ~ and is anisotropic. In the absence of 



a static field, V ~ for small distances and 

for large distances, and the interaction can be isotropic 
or anisotropic, depending on the Rydberg level. In this 
paper, we are interested in the long-range type of cou- 
pling {V ^ ^/R^)- However, to be able to simulate 
large systems, we approximate the long-range coupling 
as a constant all-to-all coupling with suitable normal- 
ization; this approximation is appropriate for a two or 
three-dimensional lattice for the system sizes used here. 

Consider a system of N atoms continuously excited by 
a laser from the ground state to a Rydberg state. Let 
\g)j and \e)j denote the ground and Rydberg states of 
atom j. The Hamiltonian in the interaction picture and 
rotating- wave approximation is {h = \) 
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where A = — Wo is the detuning between the laser and 
transition frequencies and 17 is the Rabi frequency, which 
depends on the laser intensity. 

The Rydberg state has a finite lifetime due to sponta- 
neous emission and blackbody radiation. When an atom 
spontaneously decays from the Rydberg state, it usually 
goes directly to the ground state or first to a low-lying 
state [33| ; since the low- lying states have relatively short 
lifetimes, we ignore them. In addition, blackbody radia- 
tion may transfer an atom from a Rydberg level to nearby 
levels, but this is minimized by working at cryogenic tem- 
peratures 3j|. Thus, each atom is approximated as a 



two-level system, and we account for spontaneous emis- 
sion from the Rydberg state using the linewidth 7 [24j . 
Note that each atom emits into different electromagnetic 
modes due to the large inter-particle distance. 

The environment absorbs all the spontaneously emit- 
ted photons, so the atoms are continuously monitored by 
the environment. We are interested in the temporal prop- 
erties of the emitted photons. There are two equivalent 
ways to study such an open quantum system. The first 
is the master equation, which describes how the density 
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matrix of the atoms, p, evolves in time: 
p = -i[H,p] 

+7Ef-Ji|e)(e|.'P} + l5)(eUle)(5l/) -(2) 
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A master equation of this form has a unique steady-state 
solution [35| . p^s, which can be found numerically by 
Runge-Kutta integration. The integration can be vastly 
sped up by utilizing the fact that the atoms are symmet- 
ric under interchange due to all-to-all coupling; the com- 
plexity is then 0{N^) instead of 0(4^). Using pss, one 
can calculate the statistics of the emitted light. In partic- 
ular, the correlation of photons emitted by two different 
{E^E,)/{E,){Ej), where E, = |e)(e|, [H]. 



atoms is g^^-* 



If glj' > 1, the atoms tend to emit in unison (bunching); 

('2) 

if (j^j- < 1, they avoid emitting in unison (antibunching) . 

The second approach is the method of quantum tra- 
jectories, which simulates how the wave function evolves 
in a single experiment 37-3^. In the simulation, the en- 
vironment observes at every time step whether an atom 
has emitted a photon, and the wavefunction is updated 
accordingly; the crucial point is that even when no pho- 
ton is detected, the wave function is still modified. The 
algorithm is as follows. Given the wave function \ip{t)), 
one randomly decides whether an atom emits a photon 
in the time interval [f, t + 5t\ based on its current Ry- 
dberg population. If atom j emits a photon, the wave 
function is collapsed: \'il){t + 5t)) = \g){e\j\'4'{t)) ■ If no 
atoms emit a photon, |'0(t + 5t)) = (1 — iHcS&t)\'il){t)) , 
where ifeff = H — {ij/2)J2j |e)(e|j. After normalizing 
the wave function, the process is repeated for the next 
time step. The nonunitary part of Hcs is a shortcut to 
account for the fact that the non-detection of a photon 
shifts the atoms toward the ground state 37|, |38| . 

These two approaches are related: the master equation 
describes an ensemble of many individual trajectories 37, 
HI]. Also, Pss can be viewed as the ensemble of wave 
functions that a single trajectory explores over time. We 
will use both approaches below, although quantum jumps 
are most clearly seen using quantum trajectories. 

We first consider the case of = 2 atoms since it is in- 
structive for larger N. Laser excitation and spontaneous 
emission distribute population through out the Hilbert 
space, {\gg), \ge), \eg), |ee)}. When A = 0, |ee) is uncou- 
pled from the rest of the states due to its energy shift, 
so there is little population in it [Fig. H^a)]; this is the 



well-known blockade effect [13|, [18|. But when A 



there is a resonant two-photon transition between \gg) 
and |ee), so |ee) becomes populated [Fig.[T]Jb)]. Using the 
master equation, one can calculate the photon correlation 
between the two atoms (Figs. There is strong anti- 
bunching for A and strong bunching for A « V/2, 
which makes sense since a joint emission requires pop- 
ulation in lee). One can solve for the correlation as a 
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FIG. 1: Two atoms, (a) When A = 0, |ee) is uncoupled from 
the other states, (b) When A = V/2, there is a resonant 
two-photon transition between \gg} and |ee). (c) Quantum 
trajectory simulation with Q = 1.57, ^ = = 57, showing 
Rydberg population of each atom over time. Atom 1 (solid 
blue line) emits at t — 14.4/7, which causes {E2) (dashed red 
line) to suddenly increase. Atom 2 then emits at t = 14. 7/7. 
When no photons have been emitted for a while, the wave 
function approaches a steady state. 




FIG. 2: Photon correlation for two atoms with V = IO7. (a) 
Correlation vs A for = O.57. (b) Correlation as a function 
of Q and A using color scheme on right. 
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Note that the correlation can be made arbitrarily large 
by setting w 0, A = V/2, and V large; this may be 
useful as a heralded single-photon source [40]. 

Further insight is provided by quantum trajectories. 
An example trajectory for A — V/2 is shown in Fig.[ljc). 
The atoms emit photons at various times. When no pho- 
tons have been emitted for a while, the wave function 
approaches an entangled steady state due to the balance 
of laser excitation and nonunitary decay from the non- 
detection of photons [4l| : 



cilgg) + C2\ge) + cz\eg) + Ci\ee) 
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where the coefficients have constant magnitudes and their 
phases evolve with the same frequency (this is a periodic 
steady state). Because of the laser detuning, |cip is much 



|c4| , which are comparable to each 



larger than |c2p, [csl 
other. Thus {Ei), {E2) ~ and the atoms are unlikely 
to emit. But when atom 1 happens to emit, the wave 
function becomes 
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FIG. 3: (a) Fixed points of mean-field model as function of 
detuning for Q = I.57 and V = IO7. Stable (unstable) fixed 
points are denoted by solid (dashed) lines, (b) Mean-field 
bistable region (black) for V = IO7. (c) Photon correlation 
gr'^' for 16 atoms with same parameters as (b), using color 
scheme on right. 
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Now, {E2) is large and atom 2 is likely to emit, which 
leads to photon bunching [Fig. [TJc)]. 

Then we consider the case of large N. We first review 
mean-field theory, since it is important for what follows 
24, 4^]. Mean-field theory is a classical approximation 
to the quantum model: correlations between atoms are 
ignored, and the density matrix factorizes by atom, p = 
Pj where p evolves according to 



Pee = -f^Im p^g ~ 
- - i(A - Vp^^)p^g - 
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These are the optical Bloch equations for a two-level 
atom, except the effective laser detuning is A — T^Pgg. 
There are one or two stable fixed points, depending on 
the parameters (Fig. |3]). Classically, the system should 
go to a stable fixed point and stay there, since there are 
no other attracting solutions. 

Now we consider the original quantum model for large 
N. Figure mja) shows a quantum trajectory for N — 16. 
We plot the average Rydberg population of all the atoms 
(E), where E = J^i Ei/N. (E) appears to switch in time 
between two values. In fact, these two values correspond 
to the two stable fixed points of mean-field theory for 
the chosen parameters. Thus, we find that the quan- 
tum model is able to jump between the stable states of 
the classical mean-field model. When the parameters are 
such that mean-field theory is monostable, (E) remains 
around one value and there are no jumps. As a result, the 
photons are bunched when mean-field theory is bistable 
and uncorrelated otherwise. This correspondence is ev- 
ident in Fig. [2Ib)-(c), with better agreement for larger 
N. 

We call the two states in Fig. SJa) the dark and bright 
states, since the one with lower {E) has a lower emission 
rate. In the dark state, the wave function approaches a 
steady state, \ip)ss, in between the sporadic emissions. 
This is due to the balance of laser excitation and non- 
unitary decay from the non-detection of photons, similar 
to the case of two atoms. In the bright state, the large 
Rydberg population brings the system effectively on res- 



FIG. 4: Quantum trajectory of 16 atoms showing average 
Rydberg population over time with fl = 1.67, V = IO7, A = 
3.47. (a) Quantum jumps between two metastable collective 
states. Red arrows point at the stable fixed points of mean- 
field theory, (b) and (c) are zoomed-in views, and red lines 
mark photon emissions, (b) Rapid succession of emissions 
around t — 232/7 causes a jump up. (c) Absence of emissions 
around t = 313/7 causes a jump down. 



onance. The bright state sustains itself because an atom 
is quickly reexcited after emitting a photon. 

Suppose the system is in the dark state. The steady- 
state wavefunction \ip)ss is an entangled state of all the 
atoms with most population in \gg ...g). Although {E) 
is small, when an atom happens to emit a photon, {E) 
increases due to the entangled form of |V')ss- In fact, 
if more atoms emit within a short amount of time, {E) 
increases further [Fig. [5ja)]. When enough atoms have 
emitted such that {E) is high, the system is in the bright 
state and sustains itself there [Fig.|4l[b)]. If too few atoms 
emitted, the system quickly returns to \i))ss- 

Then suppose the system is in the bright state. There 
are two ways to jump to the dark state: most of the 
atoms emit simultaneously or most of the atoms do not 
emit for a while (the non-detection of photons projects 
the atoms toward the ground state). For our parameters, 
simulations indicate that the latter is usually responsible 
for the jumps down [Fig. |4l^c)]. 

The jumps are inherently collective, since they result 
from joint emissions or joint non-emissions. As N in- 
creases, the dark and bright periods become longer and 
more distinct [Fig. [5l^b)-(c)]. This can be understood in- 
tuitively as follows. Suppose the system is in \'>p)ss- As 
N increases, the increment of {E) per emission decreases 
[Fig-E^a)]. Thus, for large iV, a rapid succession of many 
emissions is necessary to jump to the bright state. Al- 
though the emission rate in the dark state increases with 
N , the rate of nonunitary decay in Hcs also increases 
with N . The result is that the probability rate of a jump 
up decreases. Then suppose the system is in the bright 
state. As N increases, a jump down requires more atoms 
to not emit in some time interval, so the probability rate 
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FIG. 5: Statistics for SI — 1.5"f, A = 8.67, V — IO7 comparing 
N = 4,8, 16. (a) Rydberg population (E) after a number of 
simultaneous emissions from the steady-state wave function 
\tp)s3- (b) Length distribution of bright periods (in units of 
1/7), using arbitrary threshold of {E) = 0.2 and sampling 
rate of IO7. (c) Distribution of {E). Red arrows point at 
stable fixed points of mean-field theory. 



of a jump down decreases. 

These collective jumps are reminiscent of a familiar 
classical effect. It is well known that adding thermal noise 
to a bistable classical system induces transitions between 
the two stable fixed points 43|,|44|. In contrast, the jumps 
here are induced by quantum noise due to entanglement 
and quantum measurement. We note that the jumps 
may be the many-body version of quantum activation, 
in which quantum fluctuations drive transitions over a 
classical barrier 
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Experimentally, the jumps may be observed in a 2D 
optical lattice of atoms with a static electric field normal 
to the plane for long-range Rydberg interaction. For ex- 
ample, two ^''Rb atoms in the \n = 15, g = 14, m — 0) 
Rydberg state have a coupling of about 44 kHz at a dis- 
tance of 13 fim IJ], and the linewidth is 7/27r « 68 kHz 
at K [131 ■ This corresponds to V ^ IO7 with TV = 16 
atoms for the all-to-all model in Eq. (jlj, which are the 
parameters used in our discussion. One could observe 
the jumps directly by monitoring the fluorescence from 
the atoms. Alternatively, one could make repeated pro- 
jective measurements and thereby infer the existence of 
two metastable states from the distribution of (E) . 

Thus, atoms coupled through the Rydberg interaction 
exhibit collective quantum jumps. It would be interest- 
ing to see whether similar jumps appear in other settings, 
such as coupled optical cavities |25rt28i | and quantum- 
reservoir engineering 2^32|- In particular, since mean- 
field bistability seems to predict collective jumps in the 
underlying quantum model, one should look for bistabil- 
ity in the mean- field models of other systems [28, 3^, 31 1. 
Finally, we note that one can observe the conventional 
type of quantum jump in a three-level atom by using the 
Rydberg level as the metastable state [3|- Then due to 
the Rydberg interaction, a jump in one atom will en- 
hance or inhibit jumps in its neighbors. This may lead 
to interesting spatiotemporal dynamics. 
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